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Abstract 

We study integral functionals constrained to divergence- free vector fields 
in on a thin domain, under standard p-growth and coercivity assump- 
tions, 1 < p < oo. We prove that as the thickness of the domain goes to 
zero, the Gamma-limit with respect to weak convergence in is always 
given by the associated functional with convexified energy density wherever 
it is finite. Remarkably, this happens despite the fact that relaxation of 
nonconvex functionals subject to the limiting constraint can give rise to a 
nonlocal functional as illustrated in an example. 
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1 Introduction 



This article is devoted to the study the "effective" value per unit volume of func- 
tionals constrained to solenoidal (i.e., divergence- free) vector fields defined on a 
thin domain a; x (0,£), in the limit as the thickness e goes to zero. We assume 
that on a domain with finite thickness, our functional (which we call the "energy", 
although its meaning might be different from a physical point of view) is given 
by a integral of the form 

, - / 9{y',v{y)) dy if f e V^, 

Geiv) = \ ^ Jujx{0,e) 
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where > 3, a; is a bounded domain in R^"-*^, y — {y',yN) £ x (0,£), g : 
oj X ^ M is a given energy density, and is finite only in the class of 
solenoidal vector fields on a; x (0, s) in for some 1 < p < oo, i.e., 

:= {v e LP{u X (0,£);R^) jdivu^O}. 

Here and throughout the rest of this article, differential constraints as for v above 
are understood in the sense of distributions, in particular, divv = for a v E 
U'{u X (0,e);]R^) means that /^^(oe)^ ' ^'^dy — for all test functions (p e 
C^{u! X (0, £)) (smooth functions with compact support, scalar- valued) . Using 
rescaled variables given by a; = (x'.xn) = (y', e~^yiv) and u{x) — v{x',exN), 
is transformed into a functional defined on a fixed domain: 

„ . X _ J f^f(x, u{x)) dx, ii u E Ue, with f2 := a; x (0, 1), 

where f{x, •) = g{x', •) for x = {x',xn) e IR''^"^ x R, 

Ue:^{ue LP{n;M.^) \ div^u = O} 

and 

for u — {u', u^) = (m^, . . . , u^). As this does not further complicate our approach, 
we allow / to explicitly depend on xat as well below. We assume that 

/ : X ^ M is a Caratheodory function^ (f:0) 

satisfying the following structural conditions: 

(growth) |/(x,/x)|<C|^r + C, (f:l) 

(coercivity) /(a;,/x) > ^ |//|^ - C, (f:2) 

with constants C > and 1 < p < oo, for every ^ G and a.c. x G 17. 

Using the notion of F-convergence introduced by De GlORGi [10, 9], the effective 
energy of in the limit £ — >■ O"*" is expressed by the F-limit of F^. with respect to 
weak convergence in LP. For an introduction to the theory of F-convergence, the 
reader is referred to [8] and [4]. We use the notation 

r(-^weak) ~ liminf Fe{u) := inf | liminf Fs{ue) I u weakly in L^j, 

e-5>0+ 

^(-^weak) ~ linisup Fe{u) :— inf | limsup Fe{ue) \ ^ u weakly in i/j. 

£->-0+ 

■^i.e. measurable in its first and continuous in its second variable 
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Below, we omit the topology indicated in brackets as throughout this paper, this 
is always the weak topology in U*. We say that V — limF^ exists if F — liminf 
and r — limsup coincide, in which case this quantity is denoted by F — limFg. 
In particular, the use of the weak topology in U' causes a process of relaxation 
in the limit, roughly speaking because energetically favorable microstructures of 
a characteristic size converging to zero as £ — > are allowed along the sequences 
generating the effective (macroscopic) limiting energy. 

The corresponding problem of dimension reduction for functionals depending 
on gradients instead of divergence- free fields was investigated by Le Dret and 
Raoult [18, 19, 20] and stimulated a great deal of further research, including 
the study of different scalings, partially with energy densities that are realistic 
from the point of view of hyperelasticity (see [15] and the references therein), as 
well as extensions to non-flat hmiting surfaces [23, 22]. 

Recently, dimension reduction problems for Ginzburg-Landau-type functionals, 
involving a magnetic potential which is divergence-free as a choice of gauge, were 
studied in [6] and [1]. In both cases, the relevant parts of the energy density 
(apart from compact perturbations) are convex and thus no relaxation occurs 
during the limit process, avoiding the main difficulty of our problem. Relaxation 
and homogenization of functionals constrained to solenoidal matrix fields were 
treated in [29] and [2] (for related results and some physical background also see 
[16] and [27]), as well as in [14], [5] and [12] for a more general constraint of the 
form Au = 0. In this context, ^ is a linear differential operator assumed to satisfy 
MuRAt's condition of constant rank [26], and apart from the examples in [31], 
[24] and [21], very little is known if this condition is violated. In our framework, 
div^ satisfies the condition of constant rank for each e, but the associated limiting 
operator divg (divow := dp^u^ for u : Q M^) does not. From the point of view 
of the theory for ^-free fields developed in [14, 5], this means that important 
bounds for the projection operator onto divg-free fields and its complementary 
projection are not uniform in e and projecting tends to create large errors as 
£ — )■ 0"*" (cf. Remark 2.8). Hence, we can (and do) use the projection only along 
sequences that are asymptotically div^-free in a very strong sense (cf. Lemma 2.9). 

As we shall see, the divergence-free dimension reduction problem with nonconvex 
energy density exhibits some intriguing features that do not occur in the gradient 
case. In particular, it turns out that dimension reduction and direct relaxation 
in the limit setting do not yield the same result in general. While the former 
simply leads to convexification by our main theorem stated below, the latter 
may give rise to a nonlocal functional as illustrated by the example discussed in 
Proposition 3.3. 

Unless indicated otherwise, we assume throughout that 
N >2, uj C. is open and bounded, Q := a; x (0, 1) and 1 < p < oo. 
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Theorem 1.1. Suppose that (f:0)-(f:2) are satisfied. Then V — Yan.^_^Q+ (with 
respect to weak convergence in LP) exists, and it has the representation 



r — lim iu) 



F**{u) := J^r*{x,u)dx tfueUo, 

+00 ifueLP{n;R^)\Uo, 



where for each x, f**{x, •) denotes the convex envelope of f{x, •) and 
Uo:^{ue U'{n;R^) \ Onu^ = m Q} . 

It is fairly easy to see that both F — lim sup Fs{u) and F — hminf Fs{u) are fi- 
nite if and only ii u E Uq (Lemma 2.2 and Lemma 2.3), and the lower bound 
for F — liminf ^^(m) is of course a simple consequence of the weak lower semi- 
continuity of convex functionals (Proposition 2.6). However, the upper bound, 
F — limsupFe(u) < F**{u) for w e ZYq, is far more difficult than in the gradi- 
ent case. The main issue here is that a priori^ we do not know whether or not 
F — lim sup is a local integral functional. The usual trick for a proof of this 
property, based on "localizing" a sequence that weakly converges to zero by 
multiplying it with suitable smooth cut-off functions with the desired support, 
does not work in our setting, at least not in direction of the last variable, because 
the distance of the modified sequence to the set of divg-free fields in may be 
of an order approaching 1/e which is an error too large to handle. Indeed, our 
proof of the upper bound in Section 4 (culminating in Proposition 4.9) does not 
use this kind of truncation in direction xjv, instead relying on a rather explicit 
construction of suitable sequences with small support in direction of a; at which 
are asymptotically diVg-frcc in the sense that their distance to with respect 
to the norm of goes to zero as £ — )> 0+ (by Lemma 2.9). A prototype of this 
construction for a simple example is presented in Proposition 3.5. 



2 Preliminary observations 

We first observe that both F — lim sup Fs{u) and F — lim inf -Fe('u) are finite if and 
only if M e Uq. The following simple density result turns out to be useful. 

Lemma 2.1. With respect to the strong topology in i7(Q; R^), Uq n C°°(Q; R^) 
is dense in Uq. 

Proof. Let u G Uq, and extend u = {u^, . . . ,u^) to a function in LfQ^(r2;M'^) 
such that 14^' = on \ Q for j = 1, . . . , AT - 1, = on \ (w X M) and 
u^{x', xn) is still constant in x^ for a.e. x' G cu. Mollifying in the usual way yields 
a sequence {uk)keN in C°°{R^; M^) n Uq with Uk ^ u strongly in L?'(Q; M^). □ 
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Lemma 2.2. Let Vn be a bounded sequence in L^(Q) with Vn v^o weakly in 
LP{Q), and suppose that d^Vn in the sense of distributions. Then Voo is 
constant in xn- In particular, if {ue) C Ue and ^ u weakly in LP , then 
u e Uq. 

Proof. For every (p e C^{uj) and every r/ e C^((0, 1)), we have 
0= lim / / Vn{x',XN)ip{x')rj{xN) dxjsfdx' 

Jcu J {0,1) 

= / / Voo{x',XN)^{x')rj{xN)dxNdx'. 

Ju) J {0,1) 

In particular, since (p was arbitrary, we get 

/ Voo{x' ,XN)f]{xN) dxN = for a.e. x' & u and every r] e C^((0, 1)), 

J{0,1) 

which in turn imphes that Voo{x',xn) is constant in x^. □ 

Lemma 2.3. For every u & Uo, there exists a sequence (u^) C Ue such that 
Ue-u^O in ]J'{n;R^). 

Remark 2.4. Using Lebesgue's theorem, (f:0) and (f:l), we get that lim Fe{ue) — 
f{x, u) dx < oo, and thus F — hmsup Fe{u) < oo for every u eUq. 

Proof of Lemma 2.3. Step 1: Assume in addition that u e (7^(0; M^). 
For J = 1, . . . , — 1 define ui :— , and let 

/•XM 

u^{x\xn) '■= u^{x',xn) — £ / div' u{x' ,t) dt, 

Jo 

where div' u — div} + . . . + Bn-iu^'^. We thus have that divg u^ — ^ and ^ u 
strongly in LP, whence Ve '■— has the asserted properties. 

Step 2: The general case. 

By Lemma 2.1, there exists a sequence (itfe) C C^(Q;R^) fl Wo with Uk ^ u 
strongly in L^(Q; M.^) as A: — > oo. For each k and each e, we define Uk,e G ^4 as in 
the first step, using Uk instead of u. Now choose {k{e))s>o with k{e) — )■ oo slow 
enough such that e ||^fc(£) Hc-ii-q.jjjv) — )■ as e — )■ 0. As a consequence, -Ug := Uk(£),£ 
converges to u strongly in L^, and it satisfies div^ = by construction. □ 

To prove the lower bound F — hminf Fs{u) > F**{u) for u e Uq, we first recall the 
well known characterization of weak lower semicontinuity of convex functionals: 



5 



Theorem 2.5 (see [17] or [13], e.g.). Suppose that f satisfies (f:0). Then the 
functional J : U>{n,R^) [0,oo], J{u) := J^f{x,u) dx, is lower semicontin- 
uous with respect to weak convergence in LP if and only if f{x, •) is convex for 
a.e. X & fl. 

As an immediate consequence, we have 

Proposition 2.6 (lower bound). Suppose that the assumptions of Theorem 1.1 
hold. Then for every u &Uq, 

r-liminf F^(k) > F**{u). 

For the upper bound, we have to construct a suitable sequence {u^) C such 
that Ue ^ u in LP and F^i^Ue) — )■ F**{u), starting from a given u e Uo- The main 
problem here is the constraint div^ — 0. In particular, we rely on a projection 
onto divg-free fields, which is based on the following special case of the projection 
used in [14]. 

Lemma 2.7. Let 1 < p < oo and let Q C he an open cube. For every £ > 0, 
there exists a linear operator Ve '■ LP{Q]W^) — >■ U'{Q]W^) with the following 
properties: 

(i) dlY^VeU = on for every u e LP{Q]M.^), where VeU is extended Q- 
periodically. 

(a) VeW — w for every w e LP{Q; R^) such that divg w — on , where w is 
identified with its Q -periodic extension to . 

(Hi) ||^£w||^p(Q.]giV) < C'e ||w||j;,p(Q.]Rjv) for every u e LP^Q-^M.'^), with a constant 
Ce > independent of u. 

(iv) \\{L -Ve)u\\j^p^^Q.^N) < Ce\\d\Yeu\\w-'^.p{Q) for every u e Lp{Q]R^), with a 
constant Cs > independent of u. 

Here, on a given domain W~^'^ denotes the dual space ofW^'^ withp' = p/ (p— 1). 

Proof. For ^ = (^', ^^) e \ {0}, (^', ^ M^^^ has full rank independent of 
{ 7^ 0, which means that for fixed £, div^ satisfies Murat's condition of constant 
rank ([26]). Hence, Lemma 2.14 in [14] apphes with A :— div^ and T = Ve. □ 

Remark 2.8. If p = 2 (avoiding the use of general Fourier multiplier theorems), 
it is easy to see from the proof of Lemma 2.14 in [14] that (iii) and (iv) actually 
hold with a constants independent of £. However, we do not exploit this fact, 
and in any case, the factor ^ hidden in the div^ on the right hand side of (iv) is 
still a major obstacle even if the constant in (iv) does not blow up as £ —> 0"^. 
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For technical reasons, it is important for us to be able to work with sequences 
which are not div^-free but can be projected to div^-free sequences with an error 
that is negligible in the limit e — )> 0^. The following application of Lemma 2.7 
gives a useful sufficient criterion for sequences with this property. 

Lemma 2.9. Let Vt C he open and hounded, let 1 < p < oo and let En O"*". 
Then there exists a sequence an 0+ such that the following holds: For every 
sequence (un) C L^(f2;M^) with Un ^ in and 

II div.„ linL-i.,(n) + II K, ^^n ) L-i..(n;M^) < '^n, (2-1) 

where u'^ := (m^, . . . there exists a sequence (vn) C L^(r2;M^) such that 

divj^ v„ = m Q and Un — Vn ^ in LP{fl; M^) . 

Proof. For every k eN choose a function cpk G C~(Q; [0, 1]) such that (pk{x) — 1 
for every x & Q, with dist (x; dQ) > ^. Moreover, choose a cube Q containing Q 
and a sequence 5"„ — )■ O"*" such that C^^a^ with the constants of Lemma 2.7 
(iv) (which also depend on Q). We define 

C^n ||v^i(«)||vv'2,oo(f^) 5-„ and Un (Pj{n)Un 

with a sequence of integers j{n) oo (fast enough) such that m„ — ^ in 
LP(0;M^). Since 

we have that 

\\diy,^{cpkUn)\\w-i,v < lkfc|lw^2,oo(n) ( ||diVe„ 1in|U-i,p + || «, ^t^n ) IU-lp) ■ 

Hence, (2.1) implies that 

Csn l|diVe„ltn|U-i,p(n) = C'e^ ||diV£„{ln||^_i,p(Q) < Ce„a„ 

as n — 7- oo. The sequence f„ := 

VnUn e LP{Q;R^), restricted to now has the 
desired properties by Lemma 2.7. □ 

Applying Lemma 2.9 is not easy because (7„ might converge to zero extremely 
fast. Nevertheless, it turns out to be possible for certain sequences constructed 
below, first in Proposition 3.5 for a simple example and then in Proposition 4.3 
as the first step in proof of the upper bound. 
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3 An example and a related relaxation problem 



When studying the dimension reduction problem for functionals depending on 
gradients (instead of divergence- free functions), one usually relies on a character- 
ization of the associated relaxed functional in the limit setting, both as a lower 
semicontinuity result for the lower bound and as a first step in the construction 
of a sequence for the upper bound. In our framework, the associated relaxed 
functional in the limit setting corresponds to the functional Fq introduced below. 
Although Fq does not play a role in the proof our main result, we briefly discuss 
it here to point out the somewhat surprising fact that Fq does not always give 
the right limiting model for the divergence-free dimension reduction problem and 
may even be nonlocal, in sharp contrast to the gradient case. In addition, the 
crucial idea for the proof of the upper bound in our main result is developed in 
Proposition 3.5 for a simple model problem. 

In the following, we consider the functional 

F(u) ■= I ■^^ ""^ " ^ 

^ ' ' \ +00 if M ^ Uq. 

By definition, the relaxed functional associated to F is given by the lower semi- 
continuous hull of F with respect to weak convergence in L^. For u e Lp{Q; IR^), 
it can be expressed by 

Fq{u) := r — lim. F{u) — inf | lim inf F{un) Un ^ u weakly in L*'| (3.1) 

Here, note that since F does not depend on n, F — lim inf F = F — lim sup F. 
Moreover, Uq is weakly closed in L^, whence Fq{u) is finite if and only Hue Uq. 

Proposition 3.1 (partial representation of Fq). Let f : cu x — j- [0, 00) (iden- 
tified with f : D, X M.'^ — )■ M constant in xn) satisfy (f:0)- (f:2). Then for every 
u e U'{uj\M.^) (identified with u e L*'(Q;]R^) with all components independent 
of xn), we have Fq{u) — F**{u), the convexified functional. 

Proof. Since / > /** and F** is weakly lower semicontinuous in L^(r2;]R^), it 
is clear that Fq{u) > F**{u). On the other hand, for any u e LP {9,; R^) which is 
constant in xn, we have 

Fo{u) = inf jliminf F(m„) Un ^ u weakly in L^iQ; R^), OnU^ = G m| 

< inf jliminf F(k„) Un ^ u weakly in L^'(0;K^), ^jv^n = e M^} 

Un ^ u weakly in L^(a;; R^) 



inf I lim inf j f{x', Un) dx' 
I f**{x',u)dx' ^ I f**{x',u)dx, 
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where we used that f**{x' ,v) dx' is the weakly lower semicontinuous hull of 
^ ^ L /(^'' dx' in LP. □ 

Example 3.2. Let p = 6, let iV = 2, let / : M be the three-well potential 

given by 

/(/i) := l/i - Cil^ l/i ~ (2!^ 1^ - Csl^ , 
withCi := (0,-1), (2:= (1,0), (3:= (0,1), 

and consider the function uq G Uo given by 

'^'^^'^^''^^^-t (1,0) ifa;2e(i,l). 

Proposition 3.3 (Possible nonlocal character of Fq). In the situation of Exam- 
ple 3.2, we have that 

F„(«.) > ^ ^i^F„((0, 0)) + 0)). 

In particular, Fq{u) cannot he written in the form V{u) dx with some function 
V ■.R'^ -^R, and Fo(mo) > F**{uo). 

Remark 3.4. As recently discovered in [7], the lower semicontinuous hull with 
respect to strong convergence in of certain integral functional of the form 
u ^ J^f{u,'Vu)dx can also be nonlocal, if there is a lack of cocrcivity with 
respect to the gradient variable. 

Proof of Proposition 3.3. Since /** = on the closed triangle formed by ^1, ^2 
and C3, ^o((0,0)) = Fo((l,0)) = by Proposition 3.1. To prove that Fo(mo) > 0, 
we proceed indirectly. Suppose that -Fo('^o) = 0. By a standard diagonalization 
argument, we may choose a sequence m„ G Uq with Uq weakly in L^(r2, R?) 

such that Fq{uq) — \miF{un)- By passing to a subsequence (not relabeled), we 
may assume that Un generates a Young measure i/^, which for a.e. x G is a 
probability measure on M^, and by the fundamental theorem for Young measures 
(see [3], [25] or [13], e.g.), also exploiting that / > 0, we get that 

= Fo(mo) = lim / f{u^)dx > ! I f{i)dv,{i)dx. 

Since / vanishes only on {Ci, (^2, Ca}, ^^is implies that is supported in {Ci, (^2, Ca} 
for a.e. x, i.e.. 
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where 5^ denotes the Dirac mass concentrated at the point 2; in R^. Moreover, 
since /]g2 d^xiO = ^0(3;) and is a probabihty measure for a.e. x, the coeffi- 
cients <Jj{x) e [0, 1] are determined by the hnear system 

E5=i = ""0(2;) and Xl5=i (^jix) = 1- 

One easily checks that the unique solution of this system is given by 

ai{x) = i, 0-2(0;) = 0, 0-3(0;) = i iix2<l ^o(a^) = (0, 0)), 
(7i(a;) = 0, (72(x) = 1, (73(2;) = if a;2 > I (i.e., ^0(0;) = (1, 0)). 

In addition, the marginal of i/x on the second coordinate axis, 

I'liA) := ^^{u! X A) for A C (0, 1) Borel-measurable, 

is the Young measure generated by uf^ and thus independent of X2. However, 
this contradicts (3.2) and (3.3), because the latter imply that z/^ = ai{x)6-i + 
O'2{x)5o + a3{x)5i, and the coefficients are not constant in X2 (only piecewise). □ 



The dimension reduction problem is different because the constraint div^ — 
is actually genuinely less restrictive than d^u^ — 0: 

Proposition 3.5. In the situation of Example 3.2, for every given pair of se- 
quences Sn — >■ 0"^ and an — >■ 0+; there exists a bounded sequence (un) C L°°(Q;]R^) 
such that Un ^ in LP , 

I f{un + uo)dx^ I f**{uo)dx = 0, (3.4) 

and 

l|dive„ Mnllv^-i.p(n) + II K, ) IL-i,P(n;M^) < (3-5) 

for every n. In particular, Un can he projected onto Ue„ with an error that goes to 
zero strongly in W by Lemma 2.9, and since div^^ uq — div'iig = 0, this entails 
that r - liminf Fe„(Mo) < < Fq{uq). 



Proof. For each n fix a function ipn G C^((0, 1); [0, 1]) such that = 1 on 
[sn, 1 — £„], and for A; e N let 



Wk{t)^{wl{t),wl{t)) :-- 



C3 = (0,l) ifO<t<^, 
Ci = (0,-1) if ^<i<i, 



extended periodically to a function : R — > with period |. Note that 

Wk^ k3 + ki = {0,0) weakly in L^'(r;R2) (3.6) 
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for any bounded open set T C R. We define Vk,n € M^) by 

(Pn{2x2)wk{^xi) if < < ^, 



if i < X2 < 1. 



Observe that although Vk^n is not continuous, its jumps do not contribute to 
dive„ Vk,n (as a distribution), and thus the latter is actually a function with 

2 

diVenVk,n{^l,X2) = —(Pn{2X2)wl{^Xi) . 

In particular, as A; — )■ oo for fixed n, div^^ i'fc,„ ^ weakly in as a conse- 

quence of (3.6), and thus dw^_^Vk^n strongly in W~^'^{Q), by compact em- 
bedding. Analogously, we get that Vk.n — mq ^ in in PF~^'^(0;M^) as A; ^ oo. 
Hence, we may choose k — k{n) with k{n) — >■ oo as n — )■ oo fast enough such that 
(3.5) holds for Un ■— v^n),n- Again using (3.6), it is not difficult to check that 
weakly in i/, and (3.4) holds as well. □ 

Remark 3.6. The choice of the dimension N — 2 is not crucial for Example 3.2, 
it is just the simplest possible case. In fact, a completely analogous argument 
can be used for suitable potentials / with N + 1 wells in for any N > 2. 



4 The upper bound 

In this section, we provide the remaining part of the proof of Theorem 1.1, namely 
the upper bound 

r — limsup Fe{ii) < F**(u) for u e Uq, 

by constructing a suitable recovery sequence. In particular, we need some results 
from convex analysis: 

Lemma 4.1 (Caratheodory's theorem, see [30], e.g.). Let g : — )■ [0, oo) be 
continuous. Then for every ^ G and every 6 > 0, there exists an m E 
{0, . . . , A^} and e , Oj E (0, 1], j = 0, . . . , m, such that J2j = 1; C = 

and the vectors — ^q, j — 1, . . . ,m, are linearly independent. Here, g** denotes 
the convex envelope of g. 

Lemma 4.2. Suppose that the assumptions of Lemma 4-1 hold. If, in addition, 
there exist constants p > 1 and C > such that 

^ l/^r - C < gin) < C + C for every n e R^, (4.1) 
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then the assertion of Lemma 4-1 stays true even for 5 — 0, and in this case, 

<i^(|e| + l) forj = Q,...,m, (4.2) 
where K is a constant that only depends on p and C. 

Proof. With some background in convex analysis, this is not hard to prove, and 
we just sketch some details: It is well known that the convex envelope of g can 
be represented as 

g**{i) = sup { A(0 I ^ : ^ E affine and ^ < c/} , ^ e M^. 

If g is (lower semi-) continuous and has super linear growth, the supremum is 
attained at a suitable affine function (see [13], e.g.), and always touches g 
from below at suitable points as in Lemma 4.1 with 5 — Q. In addition, as a 
consequence of (4.1), we have that 

^ -C < A^{iJi) < C + C for every e co{ij] 

(the convex hull of the points j = 0, . . . , m). Clearly, the existence of an affine 
function satisfying the latter implies that co{,^j} is bounded for fixed ^, and it is 
not difficult to obtain more precise estimates that yield (4.2). □ 



The following result is the crucial step towards the upper bound for V — hm sup 
in the general case. 

Proposition 4.3. Let N > 2, let 1 < p < oo, let I C (0, 1) be an open interval 
and let En — >■ 0+. Then for every sequence Tn — > 0+ and every pair of points 

Ci) C2 £ (iTT'd numbers 71,72 £ (0, 1) such that 7^ (2 > + 72C2 = and 
7i + 72 = 1; there exists a sequence (v„) C L°°{M.^; M^) such that 

IknlLoo < max{|Ci| , IC2I}, (4.3) 

l|diVe„ ^^nlli^-i.P(a) + II {<, iJ:^^n ) L-i,P(n;M^) < (4-4) 



(4.5) 



for every n 

Vn^O in Ll^{R^] R^) as rn^ oo, 
supp(^;„) C R^-' X U,ez(sn-z + £n/'""0' 
where /'^l := {t e / 1 dist {t; dl) > e}, and 

\{vn — Cj} n U\ — > 7j \U\ \I\ for every measurable set U C M.^ (4.6) 
and j — 1, 2. 
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Remark 4.4. The assumption ^ (2 is actually obsolete. The case of equality 
is only exchidcd above because it is much simpler and will be treated separately 
in Proposition 4.6 below. 



Proof of Proposition 4-3. For each n e N fix a function ipn G C'^'l^; [0> 1]) 
such that (fn — l on /P^"! and = on M \ and define 

V'„eC°°(M;[0,l]), V'n:= J]¥'n(- + ^). 
Furthermore, for A; e N let 



Wk{t) :-- 



Ci ifO<t<7ii, 
C2 if-72|<t<0. 



extended periodically to a function : R — >■ with period |. Note that 

Wfe ^ 71C1 + 72C2 = weakly in Lf„,(M;M^). (4.7) 

With a fixed unit vector e perpendicular to Ci ~ C2) we define Vk,n £ 
L°°(R^;K^)by 

Vk,n{^) ■= ^n(^a;jv)wfc((^a;', j^xn) ■ C12) 

Observe that although x h- )■ Wk((^x', — a^Ar) •C12) i^ continuous, it is div^^-free 
(as a distribution), and thus div^^ Vk,n is actually a function with 

In particular, as — )■ cxo for fixed n, dive„ ^ weakly in L^{Q) due to 
(4.7), and thus divsn'^k,n strongly in W~^'^{Q), by compact embedding. 
Analogously, we get that Vk^n ^ in in W~^'^{Q] M^) as k ^ 00. Hence, we may 
choose k — k{n) with k{n) ^ 00 as n — >■ 00 fast enough such that (4.4) holds for 
Vn '■— Vk{n),n: S'Hd (4.3), (4.5) and (4.6) hold by construction. □ 



Caratheodory's theorem requires convex combination of up to A?^ + 1 points, but 
Proposition 4.3 only admits two points. The following elementary lemma allows 
us to handle general convex combinations by breaking them into suitable pairs 
of two. Essentially, it states that if C = ^jij is a convex combination with 
^ e if, where H is an affine hyperplane, then ^ can be rewritten as a convex 
combination of points e such that each is a convex combination of two 
of the original points, i.e., — Pij^j + Pji^f. 
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Lemma 4.5. Let m < N and let e R^, 9j e (0, 1] for j — 0, . . . ,m such that 
Yl^=o ^3 ^ 1 '^^^ vectors — $,o, j = 1, . . . ,m, are linearly independent. Then 
there exists numbers aij e [0, 1], i,j e {0, . . . , m}, such that 



= aji, = whenever Pij = 0, >, / , <^«i = 1; (4-8) 



j=0 1=0 



1=0 



ftijAi, (4.9) 



and 



m 

w/iere 



2 

ij=0 i<j 



else. 

Here, note that e [0, 1] anc/ + /^ji = 1 if {^f - - < 0- 

Proof. Let H := {y eR^ \ = i^}. Since ^ e ^ := col^i \ ] = . . . ,m} n H 
(where co A denotes the convex hull of a set A), which is a convex polyhedral 
set, ^ can be written as a convex combination of the extreme points of S. Such 
an extreme point is either given by for some j such that if — , or it is the 
intersection of H with a line segment of the form co{^j,^j}, for indices i,j such 
that ii and lie on opposite sides of H (i.e., (^f - ^^)(^f - i^) < 0). Note 
that co{^j,^i} n if = + in this case. Hence, there exist a-ij G [0, 1] 

such that aij = aji, aij = if (3ij = 0, ^j<j c^ij = 1 and (4.10) holds. Moreover, 
since aij = aji, we have that 

^ m mm 
i,j=0 j=0 i=0 

This is another way of expressing ^ as a convex combination of the points C,j- 
Since C,j — ^o, j = 1) • • • j are linearly independent, the coefficients of the 
convex combination are uniquely determined, and comparison yields (4.9). □ 

Combining multiple instances of Proposition 4.3 with Lemma 4.5, we obtain 

Proposition 4.6. Let N > 2, let 1 < p < oo, let J C (0, 1) be an open interval 
and let — >■ 0"^. Moreover, let m < N , let ^j e and 9j e (0, 1], j = 0, . . . , 
be such that 
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and the vectors — ^o, j — 1, . . . ,m, are linearly independent. Then for every 
sequence an — >■ O"*", there exist sequences (yn), (zn) C L°°(R^;]R^) such that 

IbnlLoo < maxj l^jl and < max^ , (4.11) 

9jv|/f = div'y; = onM^, (4.12) 

l|diVe„ ZnWw-i,.^^) + II (4, ^Z^) L-i.P(f2;M^) < (4.13) 

for every n 

Vn^O, Zn^O in 1/(00 X J; R^) as £ ^ 0+0, 

supp(|/„) U supp(2;„) C R"^"^ X Kji for a compact set Kn C J 

and 

\{yn + Zn — Cj} ^ U\ — > 9j \ U\ , for every measurable U C M^"^ x J (4.15) 

e->0+ 

and every j & {0, . . . , m}. 



Proof. Let a^j and be as in Lemma 4.5, and divide the unit interval (0, 1) into 
pairwise disjoint open subintervals lij, < i < j < m (some possibly empty). 



such that \L 



aij. For £ > let 



For i < j, we define bounded sequence {yij,n)n, [Zi' 

yij,n '■— XTij{e„)^ijj 



if i = j. 

C L°°(M^;M^) as follows: 



where XTij{en) denotes the characteristic function of the set Tij{6n)- For every j, 
we set Zjj^n '■— 0. For i < j, let „ be the sequence obtained in Proposition 4.3, 



(m+l)(m+2) 



applied with / := lij, Tn : 

and 72 := — 1 — j3ji. In particular. Proposition 4.3 gives that 

\{zij,n = 6 - ^ij} n c/| 

|{%,n = -4} n^l 

for every measurable U C R'^, and 







\u\ 




u\ 


n— >-oo 






1 A, 

n— >oo 




\u\ 




\u\ 



(4.16) 



supp(^, 



C 



t,N-l 



Uk&i^nk + ejlf) with a compact J^;"' C (4.17) 



for every i < j (for i — j, (4.16) and (4.17) are trivial). In addition. 



||dive„%,„| 



w-^'P{n) 



iz'- ^v^- ) 



< 



W-^'P(n; 



(m+l)(m+2) 



(4.18) 
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for every n and every i < j. Now let 

m j m j 

z„{x) := ^^Zij^n{x) and yn{x) := ^^yij,n{^) for x e 

j=0 i=0 j=0 i=0 



Note that at any given x, at most one term contributes in each of the double 
sums above; more precisely, z^i — z^j^^ and y„ = yij^n on ^ij(£n)- Moreover, 



m 3 

£->0+ 



j=0 i=0 i<j j 



weakly in i7(Q;R^), and 

d^y^ = div'y'^ = onR^ 

since xn) is constant for every xat £ R and = = a.e.. By (4.18), we 
obtain that 



II dlV£„ 5n||^-i,p(f^) + II {Zn, er,^n) |liy-i,P(n;]RiV) 

— Z^i=0 Z^i=0 (m+l)(m+2) '^^^ ^ "^^"^ 2 ^ 

for n e N. By (4.16), we get that 

\{yn + Zn = ij) ^ U\ + ^ (3ijai^ \U\ = 9j \U\ (4 20) 

for every j and every measurable U C M^, where the latter equality is due to 
(4.9) combined with the fact that = ajj = PjjCtjj. Finally, define 

Zn ■= XR'^-'^xK„Zn and yn := XR^-^xKnVn 

where 

Kn := U {snk + e„[0, 1]) and Z„( J) := { A; e Z | £„A; + £„[0, 1] C J} . 

fcGZ„(J) 

Clearly, (4.11), (4.14) and (4.15) are satisfied, the latter as a consequence of 
(4.20). In addition, 

Zn^O and = = in a vicinity of K^"^ x dKn, 

the former by (4.17). Consequently, dnyn = ^Nlln — ^^d div'i/^ = div'^^ = 
on R^, and (4.19) implies (4.13). □ 

The next result essentially yields the upper bound in the piecewise constant case. 
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Proposition 4.7. Let be a function satisfying (f:0)-(f:2) and let & Uo- 

Moreover, let J^. C (0,1) be a finite number of pairwise disjoint open intervals 
covering (0, 1) up to a set of measure zero, let Uh G u be a finite number of open, 
pairwise disjoint sets covering cu up to a set of measure zero, and suppose that 
for each {h, k) and each n e R^, 

and f^{-,iJ,) are constant on Qh,k, where Qh,k '■— i^h x Jk- 

Then for every pair of sequences Sn — )■ 0"^ and — )■ 0"*"; there exist two sequences 
(vn), (wn) C L°°(n;K^) such that v„ ^ and w„ ^ m i;'(Q;M^), 

\vn{x)\ < K{\u#{x)\ + l) and \wn{x)\ < K{\u#{x)\ + l) fora.e.xeQ, (4.21) 

where K is a constant that only depends on the constants in (f:l) and (f:2), 

diVe^Vn^O onW^, (4.22) 

l|dive„ Wn\\w-i,P(n) + II 7:^n) Ww-^-nm^} - ^'^ ("^-^^^ 
for every n e and 

lim / fjt{x, u + Vn + Wn) dx = f#{x, u) dx, (4.24) 

where for every x, f^{x, •) denotes the convex envelope of f^{x, •). 

Proof. Step 1: We first show the assertion with (4.22) replaced by the condition 

9iv< = onR^ and ||div'<||^^(^) < (£„)-i (4.25) 

Clearly, it is enough to define Vn and Wn on each Qh,k and prove the asserted 
properties with Qh^k instead of Q, as long as the restriction of Vn and Wn to any 
one Qh,k has compact support in this set. Hence, we consider h and k to be fixed 
below. 

Let (an) C (0, oo) be a sequence with (7„ — >■ 0+ (fast enough, as specified later), 
and define 



Uh,k u#{x) and gh,k{lJ') — f#{x, A* + Uh,k) for x e Qh,k and fj, e 



By Lemma 4.1 and Lemma 4.2, G M can be written as a convex combination 
= Ylf=o ^3^3 that — C,o, j = 1, . . . ,m, are linearly independent and 

E7=o^.-/#(^, + u#{x)) = E7=o^.-^M(e,-) = 91*0) = /;*(^, u#{x)), (4.26) 
for every x e Qh,k- Moreover, as a consequence of (4.2), 

max. le,! </^(|«/^,fc| + l), (4.27) 
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with a constant K only depending on the constants in (f:l) and (f:2). Proposi- 
tion 4.6 appUed with J = Jk yields two sequences (zn) C L°°(M-'^; R-^) such 
that Hn ^ and 2;„ ^ in Lf^^, 

\yn(x)\<K{\uh,k\ + l) and \zn(x)\ < K{\uh,k\ + l) for x e M^, (4.28) 

d^y^ ^dw'y'^ = onR^, (4.29) 

l|div.„ ^n|lty-i,P(n) + II ^^n) L-i.P(f2;R^) < ^n, (4.30) 

yn and Zn vanish in a vicinity of x dJ^ (depending on n), (4-31) 

and 

» m 

lim / ghAVn + Zn) dx = \Qh,k\ Ojgh,k{Q, (4.32) 

''Qh,k j=o 

the latter due to (4.15) and Lebesgue's theorem. Together with (4.26), (4.32) 
yields that 

Jirn / gh,k{yn + Zn) dx = g**{0)dx. (4.33) 

" °°jQh,k -'Qh,k 

To obtain functions with compact support in Qh,ki we have to cut off yn and 
Zn near {duh) x Jk- For this purpose choose a sequence of functions r]n € 
C^iujh, [0, 1]) in such a way that 

r]n /^l pointwise and HVr^nll^oo < {Sn) 



K{\uh,k\ + l) 

Below, we identify r]n with a function in C°°(M^) that is constant in xn- In 
particular, we have that 

(1 - rjn)yn and (1 - rin)zn pointwise a.e. on Qh,k- (4.34) 

We define 

Vn Tlnyn and Wn := r]nZn 

By construction, these functions have compact support in Qh,ki ^ in U' and 
«;„ ^ in L^, and (4.28) entails (4.21). In addition, we have (4.25), its second 
part since by (4.29), div'v^ = (V'?7„) • y'^ = (V?7„) • yn and thus 

lldiv'v^ll^oo < llVr^nll^oo Wynhoo < (£n)~^- 

By Lebesgue's theorem, (4.33) and (4.34) yield (4.24) for Q^^k instead of fl. 
Finally, 

diYeAVnZn) = C^'Vn) ' z'^ + Vn diVe„ Zn: 



18 



whence 

^ lkn|liy2,oo(]RJV) ^ ||diV£„ ^n|lH^-i,P(n) + II {z'n, ^^)^) II vy-i,P(n;M^)) 
— ll^'re|lvF2,oc(]]j]V') C"n 

by (4.30). With := T„(||r7„||^^2,oo)~\ this gives (4.23) for Qh^k instead of fl. 

Step 2: We still have to modify f„ to obtain (4.22) instead of (4.25), while 
maintaining the other asserted properties. For x G let 

Vn{x) := Vn{x) - En div' ^^(a;', t) dt 
Jo 

for X = {x',xi\f) G Q, with Vn as in the first step. Since d^Vn = 0, we have 
diy^Vn = on O by construction, and due to the second part of (4.25), 

\\Vn - VuWl^ ^ i^n)^ 0. 

As a consequence of the latter, (4.21), (4.23) and (4.24) also hold for Vn instead 
of Vn (in case of (4.21) with a slightly larger constant). □ 

The proof of the upper bound in the general framework relies on approximation 
and the following well-known property of Caratheodory functions. 

Proposition 4.8 (ScorzarDragoni, e.g. see [11]). Let fl C M'^ be open and 
bounded and let f : x — >■ M 6e a Caratheodory function. Then for ev- 
ery 5 > 0, there exists a compact set Q C such that < S and f is 
continuous on fl x M.^ . 

Proposition 4.9 (upper bound) . Assume (f:0)-(f:2), letu G Uq and let En — >■ 0+. 

Then for every S > 0, there exists a sequence (k„) C Ue^ such that Un ^ u in 
LP {n;R^), and 

lim / f{x,Un)dx< I f**{x,u)dx + 5. (4.35) 
Jn Jn 

Remark 4.10. Since (f:2) yields a bound on ||tin|lLP independent of 5, a diagonal- 

ization argument similar to the one in the third step of the proof below shows 
that the assertion of Prioposition 4.9 stays true even for 6 = 0. 

Proof of Proposition 4-9. Using a series of approximations, the assertion is 
reduced to Proposition 4.6. Any expression of the form "A fa B" below means 
that A = B + e, with an error e whose modulus is controlled by a suitable fraction 
of (5. 
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Step 1: Assume that it e Wo is continuous in Q and / is continuous on Q x 
for some compact Q C Cl. We claim that in this case, there exists sequences 
(un) C Ue„ and (r^) C R^) such that Un ^ and r„ in U", 



hm / f(x,Un — rn)dx= hm / f{x,Un)dx!^ / f**{x,u)dx (4.36) 

and 

\un{x) - r„(x)| < {2K + l){\u{x)\ + 1) for a.e. x e Q, (4.37) 

where K is the constant in (4.28) (which, unhke Un and r„, is independent of Q). 

For the proof, we divide ^2 into sets of the form Q^^^ = '^h ^ Jk and define 
associated piecewise constant approximations of u and / as follows: Let — 
(li^, . . . ,ii#) be given by 



u 



.^(x) := inf {min{M-^(a;),0} | x G Qh,k] + sup { maxjii^ (x) , 0} | x G Qh,k] , 



for J = 1, . . . , A^, whence is a piecewise constant function in Uq such that 
< I'^'^l- Moreover, for x e Qh,k let 

■) := ■) with a fixed e { '^l^^^-l > «• 

Note that Xh,k can always be chosen in such a way that satisfies (f:l) and (f:2) 
with the original constants. In the following, let 



S := Br{0) C R^, with := (2K + l)(\\u\\^^^^.^^^ + 1), 

where K is the constant in (4.28). If the mesh size (the maximal side length of 
the boxes Qh,k) is small enough, we get that 

max\u{x) — u^{x)\ and rnax \f{x,ii)—f#{x,ii)\^0 (4.38) 

by the uniform continuity of u and / on compact sets. With the sequences Vn 
and Wn of Proposition 4.6, using (4.38), (4.21) and the uniform continuity of / 
onQ X S, we thus have that 

/ f{x,U + Vn + Wn)dx!i:i / f {x, U# + Vn + Wn) dx ^ / f^{x,U^+Vn + Wn)dx 

Jn Jn Jn 

uniformly in n. Similarly, (4.38) and the uniform continuity of /** on fix S yield 
that 

/ f**{x,u)dx^ / f**{x,u^)dx^ / f#{x,u^)dx. 
Jn Jn Jn 
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Together with (4.24) (for and instead of / and u), this gives 



n— >oo 



hm / f{x,u-\-Vn-\-Wn)dx Pi / f**{x,u)dx. (4.39) 



Finally, by Lemma 2.3 and Lemma 2.9 applied to u and Wn, respectively, there 
exists a sequence (r„) C Lp{Q.; M^) such that r„ — > in and dive„{u + Wn + 
Tn) = on Q. By Lebesgue's theorem, (f:0) and (f:l), we have that 

lim / f{x,u-\-Vn-\-Wn)dx^ lim / /(x, -u + v„ + w„ + r^) dx, (4.40) 

also using that m G is fixed and {vn)-, {w^) are bounded in L°° . Combining 
(4.39) and (4.40), we infer (4.36) for Un ■— u + Vn + Wn + r„, and (4.37) is a 
consequence of (4.28) and the fact that \u^\ < \u\ a.e. in Q. 

Step 2: Assume that u eUq is continuous in Q. 

As a consequence of Proposition 4.8, there exists a compact subset Q of Q such 
that / is continuous on Q x R^, and |n \ 0| is small enough such that 

/ \f**{x,u)\ dx^Q (4.41) 
Jn\n 

and 

sup/ \f{x,v{x))\ dx^O, (4.42) 

where V := {v e L'P{Vt] M^) | \v\ < {2K + l){\u\ + 1) a.e. }. Here, note that the 
set {fi-,v{-)) \ V G V} C L^{^1) is equiintcgrable by (f:l). With the sequences 
{un) C Ue^ and (r^) C LP{fl;M.^) of Step 1, we thus have that 



lim / f{x,u„-r„) ^ lim / f{x,Un)^ / /* 
in in in 



{x,u) 



due to (4.36), (4.41), (4.37) and (4.42). 
Step 3: The general case. 

By Lemma 2.1, there exists a sequence (uk) G Uq (1 C{Q]M.^) with ^ u in 
LP{n). Let {uk,n) C Ue^ and (ffe,„) C L^iyt- M^) denote the sequences correspond- 
ing to Uk obtained in the previous step. By (4.37), Uk,n — fk,n is bounded in L^, 
uniformly in k and n. Since the dual of 1^(0; M^) is separable, Uk,n — fk,n Uk 
in LP as n ^ oo, Uk ^ u in LP as k ^ oo, fk,n — >■ in as n — >■ oo, and 

lim / f{x,Uk,„)^ / f**{x,Uk) — > / f**{x,u), 
in in in 

there exist diagonal sequences 



Ur 



Uk{n),n e and r„ := fk{n),n e 1/(0; M^) 
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with k{n) — >■ oo slow enough such that — ^ w in I^, rn — >■ in i/, and 



Remark 4.11. It is natural to ask whether our result also holds for functionals 
on Div-free matrix fields (i.e., each column is divergence-free). The approach 
presented here extends in a straightforward way to fields with values in IR-'^^-'^ for 
M < N—1, but it does not work for M > N. Of course, for M > N, the matrices 
can have rank A^, and in general, it is no longer clear if Div-quasiconvexity (S- 
quasi convexity in the terminogy of [28], which implies convexity along directions 
of rank < — 1) implies convexity. We expect that in this case, the convex 
envelope in Theorem 1.1 has to be replaced by a suitable variant of a quasiconvex 
envelope. We hope to address this in a future work. 
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